arXiv:1504.03736vl [cond-mat.soft] 14 Apr 2015 


Margination regimes and drainage transition in confined multicomponent suspensions 

Rafael G. Henriquez Rivera, Kushal Sinha, and Michael D. Grahand 
Department of Chemical and Biological Engineering 
University of Wisconsin-Madison, Madison, WI 53706-1691 
(Dated: April 16, 2015) 

A mechanistic theory is developed to describe segregation in confined multicomponent suspensions 
such as blood. It incorporates the two key phenomena arising in these systems at low Reynolds num¬ 
ber: hydrodynamic pair collisions and wall-induced migration. In simple shear flow, several regimes 
of segregation arise, depending on the value of a “margination parameter” M. Most importantly, 
there is a critical value of M below which a sharp “drainage transition” occurs: one component is 
completely depleted from the bulk flow to the vicinity of the walls. Direct simulations also exhibit 
this transition as the size or flexibility ratio of the components changes. 


Introduction. Flow-induced segregation is ubiquitous 
in multicomponent suspensions and granular materials, 
including systems as disparate as hard macroscopic parti¬ 
cles in air (l| , polydisperse droplet suspensions (2| , foams 
3 , and blood. During blood flow, the focus of the present 
work, both the leukocytes and platelets segregate near 
the vessel walls, a phenomenon known as margination, 
while the red blood cells (RBGs) tend to be depleted in 
the near-wall region, forming a so-called cell-free or deple¬ 
tion layer [^. Engineering the margination process has 
been proposed for microfluidic cell separations in blood 
{e.g. [5l]) as well as for enhanced drug delivery to the 
vasculature 

Direct simulations of flowing multicomponent suspen¬ 
sions - models of blood - can capture margination 
phenomena [Ml, but developing a fundamental un¬ 
derstanding of underlying mechanisms and parameter- 
dependence from simulations is difficult. It is thus im¬ 
portant to have a simple yet mechanistic mathematical 
model, ideally one with closed form solutions that reveal 
parameter-dependence, that can distill out the essential 
phenomena that drive segregation and capture the key 
effects and transitions. We present such a model here. 


Theory. We consider a dilute suspension containing 
Ns types of deformable particles with total volume frac¬ 
tion (j) undergoing flow in a slit bounded by no-slip walls 
at y = 0 and y = 2H and unbounded in x and z. Quan¬ 
tities referring to a specific component a in the mixture 
will have subscript a: for example Ua is the number den¬ 
sity of component a. We consider here only simple shear 
(plane Gouette) flow and, consistent with the diluteness 
assumption, take the shear rate 7 to be independent of 
the local number densities and thus independent of po¬ 
sition. In a dilute suspension of particles, where (j) 1, 

the particle-particle interactions can be treated as a se¬ 
quence of uncorrelated pair collisions 17 For the 
moment, we neglect molecular diffusion of the particles. 
This issue is further addressed below. Since the particles 
are deformable, they migrate away from the wall during 
flow with velocity Vam{y) [13, 211. The evolution of the 


particle number density distributions can be idealized by 
a kinetic master equation that captures the migration 


and collision effects ( 3, 3 22 , 11]) • Assuming uniform 
particle distributions in x and z, this equation is 
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Here Sy and Sz are the pre-collision pair offsets in the y 
and z directions and Ay^{Sy,6z) and A'^^{6y,Sz) are the 
cross-stream and cross-vorticity direction displacements 
of a particle of type a after collision with a particle of 
type /3. See Fig.[T]for a schematic. The term 7|i5y| in the 
integrand accounts for the relative velocity of approach 
of two colliding particles. 
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FIG. 1. (a) Suspension of particles in a slit under simple shear 
flow, (b) Pair collision trajectories of particles of species a and 
/3 under simple shear flow. 


We now construct an approximation to this model 
that is valid in the limit (A“^,A“^) 0. Taylor- 

expanding the first term in the curly brackets in Eq. [T] 
about A“^ = A“^ = 0, neglecting terms involving 

(A«/3)3 smaller, and applying the condition that n„ 
is independent of z yields a set of nonlocal drift-diffusion 
equations: 
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Here Vac is the collisional drift velocity of component a, 
while Da is its short time self-diffusivity. In the impor¬ 
tant special case of a binary suspension composed of a 
“primary” component (a =‘p’) and a “trace” component 
(a =‘t’) such that Up 2 > nt, only the primary component 
contributes to these quantities: 



np{y - Sy)AyP{Sy)'^\Sy\ dSy, 


(3) 


The parameter Kam depends linearly on the yy- 
component of the stresslet generated by the deformable 
particle 2^, 211. This scales as where Oq. is the par¬ 
ticle radius of species a, and as 7 ^ at low 7 with this 


dependency becoming weaker as 7 increases [16 

With these further idealizations, Eq. [2] becomes a pair 
of partial differential equations, which we present here in 
nondimensional form: 


Da{y) = 
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outnumber platelets and white blood cells by one and 


three orders of magnitude respectively 24 1 . 


Finally, a further simplification allows substantial ad¬ 
ditional insight. We make local approximations to the 
integrals, Eqs.|3]and|H based on the argument that 
and APP are vanishingly small for large |(5y| by Taylor- 
expanding Up around 5y = 0, noting that Ay(6y,5z) is 
odd in 6y and keeping only the leading terms. Now the 
collisional drift velocities and diffusivities become 


= -K„ 


djup 

dy 


Da = K, 


oidTi'^P') 


(7) 


9(/)p 


(4) 

dt dy 

f^pm 

in- 

-K ^6 - 

Qy dp 

(5) 

d^___^ 

dt dy 

^tm j 

( 6 ) 

d(t>p, 

- Ktc-^(Pt - 

dy 


Here (j)p = UpVp and (pt = 


1 


1 

? ■ (2C - yf 
d(t>l 

dy 

1 1 

? " (2C - y)2 
d{(j)p(t)t) 
dy 


( 11 ) 


( 12 ) 


of the primary and trace components, where Va is the 
volume per particle of component a, C = H/ap is the 
confinement ratio, i^pm — — Vpa^ ’ ^pd — ’ 

Ktm = Ktc = and Ktd = Time t is nondi- 

mensionalized with 7 “^ and y with Op. For simplicity, 
we keep the symbols t and y for their nondimensional- 
ized forms. For a single-component suspension of rigid 
particles {Ng = 1 , Kam = 0 ) a model of similar form was 
proposed by based on phenomenological arguments 
first proposed by [27j |. 
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The convergence of these integrals deserves mention. In 
the far field each particle appears as a force dipole, so 
in an unbounded domain the collisional displacements 
AyP would decay as Thus convergence of Eq. [Hlis 
unproblematic irrespective of rcut- For convergence of 
Eq. [51 Tcut must be bounded. An explicit bound is the 
slit width 2H. Furthermore, at any finite concentration 
the spacing between particles scales as . A given 

particle will effectively only collide with other particles 
within this range, while particles outside this range would 
be more strongly affected by their nearer neighbors. 

To describe the wall-induced hydrodynamic migration 
velocity, we superpose the point-force-dipole approxima- 
tions corresponding to each of the two walls 2 ^ : 
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FIG. 2. Steady state volume fraction profiles of 4>p/4>p (black 
solid line) and (f>t/4>t for various values of M. (The curves 
coincide when M = 1.) Here (j>p = 0.12, 4>pc = 0.23, C = 5.08, 
Kpm = 0.11, Kpc = 0.02, and Kpd = 0.07, resulting in Id = 1.6 
(extracted from simulation results in [l 6 l|l. For simplicity, 
K-td = Kpd and Ktc = Kpc. We vary M by changing Ktm- 


Results. An important feature of Eqs. [TT] and [ 12 ] is 
that steady state solutions with no-flux boundary condi¬ 
tions at the wall {y = 0) and centerline {y = C), can be 
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found analytically. For (j)p we find 
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where (j)pc is the volume fraction of the primary com¬ 
ponent at the centerline, rir, = —i a-nd Id is the 
nondimensional cell-free layer thickness: 


we call this phenomenon the drainage transition, since for 
M < — 1 all the trace component is completely drained 
from the bulk. If the trace component does not migrate 
(as in the case of rigid particles), then utm = 0 and 
M = —{l + Ktc/ntd)-, which is always less than —1. (This 
case is degenerate in the absence of Brownian diffusion, 
because at steady state (j)t can take on arbitrary values 
when y < Id-) 
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The black solid line in Fig. [5] shows (pp normalized by 
its mean volume fraction pp. In the unconfined limit 
C —)• oo. Id — rjp/ippcj confirming the dependence 
found earlier in scaling analyses 23|, [2^, [2g . More gener¬ 
ally, Eq. [14] analytically captures the dependence of the 
cell-free layer thickness on the volume fraction, degree of 
confinement and particle properties. 

For the trace component, the steady state solution is 
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where ppiy) is the steady state solution found above, ptc 
is the centerline volume fraction of the trace component 
and 


l^tTn ^tc T ^td 

^pm ripQ -j- O-Kpd 

Remarkably, this single quantity, which we call the 
margination parameter, determines the qualitative na¬ 
ture of the concentration profile. 

The sign of M is determined by the competition be¬ 
tween the ratio of the migration velocities of the two 
components, and the ratio of the collisional terms, 
Kta+Ktd ^ Depending on M, several distinct regimes of 

f^pc~r^l^pd 

behavior can be identified: 

(1) M > 1: the trace component is displaced further from 
the wall than the primary component: it demarginates. 

(2) 0 < M < 1: the relative concentration of the trace 
component is higher near the wall than the primary com¬ 
ponent but does not display a peak: it weakly marginates. 

(3) — 1 < M < 0: the trace component displays a peak at 
y = Id, corresponding to an integrably singular concen¬ 
tration profile: it moderately marginates. 

(4) M < — 1: here Eq. [T5| displays a nonintegrable singu¬ 
larity at y = Id- This steady state is physically unrealiz¬ 
able as it corresponds to an infinite amount of material in 
a finite region. In this regime collisional transport over¬ 
whelms migration, and the trace component accumulates 
indefinitely at y = Id, indicating strong margination. 

The black solid line in Fig. [3] shows the ratio between 
the centerline concentration ptc and the average concen¬ 
tration (pt vs. M. This falls sharply to zero at M = — 1; 
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FIG. 3. Centerline volume fraction of the trace component 
ptc scaled with the average trace concentration pt vs. M for 
varying ©. Other parameters are the same as in Fig. [2l 


For particles the size of blood cells (> Ipm) at 
shear rates characteristic of the microcirculation (10^ — 
10^s“^), Brownian diffusion is unimportant. For smaller 
particles, however, such as might be used for drug de¬ 
livery, this may no longer be true. We consider the im¬ 
pact of Brownian diffusion on trace component trans¬ 
port by adding an appropriately nondimensionalized dif¬ 
fusion term tD d’^pt/dy'^ to Eq. [121 Here D = Ds/app, 
where Db is the Brownian diffusivity of the trace compo¬ 
nent. Using typical values for blood (op = 3.9 x 10“® m, 
7 = 500 s“^) and the Stokes-Einstein relation, varying 
D from 10“^ to 10“^ corresponds to varying at from 
- 10“® m to - 10-®m. 

The steady solution for the trace component with 
molecular diffusion is 
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Molecular diffusion results in a spreading of pt to include 
the region y < Id and also renders the steady solution 
for M < — 1 integrable. It also smears out the drainage 
transition as shown in Fig. [3| 

Now consider the rigid trace particle case Ktm = 0. 
Fig. [4] shows how the steady state profile of pt varies 
with ®: Margination is weakened by diffusion. Fig. [4] 
also shows the transient evolution of pt tor tD = \ x 10“^ 
from a uniform initial condition as determined from a 
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FIG. 4. Steady state volume fraction profiles (f>t{y)/4>t for 
Ktm = 0 and various ®. Other parameters are the same as 
in Fig. 2. Also, the solid lines show the transient evolution 
for ® = 1 X 10“^: t = 0 (lightest gray), 30, 90, 270, and at 
steady state (black). 


numerical simulation using a conservative finite volume 
method. 



t = 2210, F = 0.4 



F or S 


FIG. 5. Time-averaged volume fraction of the trace com¬ 
ponents at the centerline as a function of flexibility ratio 
F or size ratio S from direct hydrodynamic simulations, 
where C = 5.08 and Id = 1.6. In case (a) <(> = 0.12 
and in case (b) 0 = {0.096,0.098,0.101,0.108,0.12} for 

S = 10.2,0.4,0.6,0.8,1.0}. 

Comparison with direct simulations. To evaluate the 
prediction of a drainage transition, we performed di¬ 
rect simulations of binary suspensions of fluid-filled non- 
Brownian elastic capsules at low Reynolds number using 
a boundary integral method (cf. [10|, llfl l29| ). Two cases 
were considered: segregation by (a) deformability and 
(b) size. The particles are all spherical at rest. Particle 
deformability is characterized by the capillary number 


Gao, = p'jaa/Ga, where p is the fluid viscosity and Ga 
is the membrane shear modulus of component a. In case 
(a) the primary component comprises 80% of the parti¬ 
cles and has Ca^ = 0.5; the trace component is stiffer, 
and we define a flexibility ratio F = GpjGt = Cat/Cap. 
The primary component in case (b) is the same as in case 
(a), but now the trace component is smaller as defined 
by the size ratio S = at/op. In this case Cat = Cap. 

Fig. [5] shows the steady-state value of (j)tc/^t as F or 
S changes. It is very similar to Fig. |3l clearly indicat¬ 
ing that the drainage transition predicted by theory is 
found in the simulations. Coincidentally, the transition 
is in the same range 0.4 — 0.6 for both S and F under the 
conditions chosen. Considering case (b) first, the migra¬ 
tion parameter Ktm scales as at constant Ca*, so as S 
decreases so does M; recall that M < — 1 for Ktm = 0. 
With regard to case (a), Ktm also decreases with decreas¬ 
ing F, and additionally the collisional displacements and 
thus Ktc and Ktd increase 0- Therefore, decreasing F 
also corresponds to decreasing M, resulting in a drainage 
transition. 



S 


FIG. 6. M vs. S for Ktm ~ (black solid line) and Ktm ~ S® 
(red dashed line); the latter corresponds to the simulations 
shown in Fig. [5l Other parameters are the same as in Fig. 2. 

Returning to the theory, using the values in the caption 
of Fig. [2] and Eq. [TFl we can determine the value of S 
corresponding to the drainage transition by finding M 
as S is varied. This result is shown in Fig. [SI where 
we use Ktm ^ (black solid line) and Ktm ~ (red 
dashed line) to represent the cases of varying at while 
keeping Gt and Ca* constant, respectively. The values of 
S corresponding to the drainage transition are 0.61 and 
0.52, respectively. The latter case corresponds to case 
(b) above, and we see that the theory result agrees well 
with the direct simulation result in Fig. [5] 

For reference to blood, the values of S for leukocytes 
and platelets with respect to RBCs are 0.9 — 2 [s^ and 
0.25 [m respectively, while F is of the order 10“^ jsij 
and 10“^ [33( respectively. Thus, case (a) here is more 
closely related to the RBC-leukocyte segregation, where 
the size ratio is close to unity, and case (b) more nearly 
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represents the RBC-platelet case, where the sizes are very 
different. From the present results it appears that both 
leukocytes and platelets would satisfy the conditions for 
drainage in simple shear. 

Conclusions. A mechanistic theory has been devel¬ 
oped for the first time to describe flow-induced segre¬ 
gation phenomena in multicomponent suspensions such 
as blood. Experimental and computational observations 
of margination and demargination are captured qualita¬ 
tively in simple closed form solutions. Several different 
margination regimes arise and a sharp drainage transi¬ 
tion is identified beyond which the trace component of 
the suspension partitions completely to the edge of the 
cell-free layer. Direct simulations corroborate this pre¬ 
diction. 

The framework presented here can be extended in 
many directions. For example, pressure-driven flow, 
which is common in microfiuidic and circulatory appli¬ 
cations, can be studied. The model can be expanded to 
include many other phenomena, including red blood cell 
aggregation and platelet adhesion. Most importantly, the 
mechanistic nature of the theory leads to substantial and 
systematic insight into the origins of margination; this 
will complement detailed simulations and experiments in 
guiding development of technologies involving blood and 
other multicomponent suspensions at small scales. 
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